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Anomalous energetics and dynamics of moving vortices 
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Motivated by the general problem of moving topological defects in an otherwise ordered state and 
specifically, by the anomalous dynamics observed in vortex-antivortex annihilation and coarsening 
experiments in freely-suspended smectic-C films[T], I study the deformation, energetics and dynamics 
of moving vortices in an overdamped xy-model and show that their properties are significantly and 
qualitatively modified by the motion. 

PACS numbers: 


Introduction: Topological defects play a central role 
in phase transitions, relaxation of generalized strain in 
ordered states (e.g., current dissipation in a superfluid, 
strain relaxation in a crystalline solid, etc.) [2], coarsen¬ 
ing dynamics after a quench into an ordered state[3], 
and appear in a broad range of physical realizations 
from superfluids and liquid crystals [1] to early universe 
baryogenesis[S]. 

Many physical systems involve topological defects mov¬ 
ing (stochastically or deterministically) through an oth¬ 
erwise ordered medium. Although it is usually tacitly 
assumed that defect’s properties (texture structure, in¬ 
teraction, dynamics, etc) are not modified by its motion, 
with the center rp simply boosted rg —>■ vo{t) by the mo¬ 
tion, there is no a priori reason for this to be the case. 
Instead, not unlike a relativistic charged particle, a mov¬ 
ing defect is defined and governed by the dynamics of 
the associated vector field, requiring a nontrivial analysis 
that is the subject of this Letter. 

Stimulated by this general question, and by the 
anomalous vortex-antivortex annihilation and coarsening 
dynamics [3] observed in freely-suspended smectic-C films 
experiments [TJ [51 [7], I explored the nature of moving vor¬ 
tices in an overdamped two-dimensional (2D) xy-model, 
applicable to a broad range of soft matter systems. In 
this Letter I report the results of these studies, that, 
with some modifications may also extend to vortices in 
a nonzero-temperature superfluid and superconductor in 
the presence of a background supercurrent or dislocations 
in a strained crystal. 

Results: Before turning to the analysis, I summarize 
the results of this study. I find that a vortex imposed to 
move with a constant velocity v in an ordered medium of 
stiffness K and damping 7, beyond a length scale 

iv = — = D/v = k~^ (1) 

7U 

exhibits a nontrivial longitudinal distortion of its stan¬ 
dard, purely transverse form|3], latter retained on length 
scale below In the steady state the resulting deformed 
vortex exhibits a parabolic comet-like tail, extending 
across the system to which most of the 27r phase winding 
is confined (see Figs{g|). While a motion-induced dis- 
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FIG. 1: A vector field corresponding to a phase 6{r) for a 
vortex-antivortex pair co-moving to the right with velocity v. 
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FIG. 2: A vortex-antivortex counterpropagating pair relevant 
to the annihilation problem. The motion-induced comet-like 
tails (that lead to a linearly diverging vortex elastic energy) 
and a suppression of the deformation between vortices (that 
leads to a weakened interaction) are clearly visible. 


tortion is not surprising, the qualitative long-scale nature 
of its consequences (see below) indeed is. For a transient 
state at time t after a vortex begins to move, the steady- 
state distortion only extends out to a time-dependent 
anisotropic “horizon” vt x ^J~KtJ^, beyond which the 
purely transverse vortex field is nearly undistorted by the 
motion. This is analogous to the Lienard-Wiechert po¬ 
tential of a moving point charge |9]. All other predictions 
follow from this result. Specifically, the vortex steady- 




2 


X \ \ t 
'v \ \ 't 
X X \ 

X \ \ 
^ X X \ 


X \ 


//XXX 
/ / X X X 
/ X X X X 
/ X X X X 
/ X X X X 
f / X X ^ 

X \ X ^ 


X / 


X 


XXX/ 
XX// 
XX// 
X / / / 
^ / / / 


/ X — 

1 \ X X 
\ X X X 
\ X X X 
\ X X X 
\ \ X X 
\ \ X X 


FIG. 3: A transient vector field of a vortex moving for time t, 
exhibiting a steady-state distortion out to an elliptical “hori¬ 
zon” ft X and purely transverse vortex field beyond. 
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vanishes logarithmically with vanishing velocity, cut¬ 
ting off the InL/a divergence of a stationary vortex drag 
coefficient, a result that was previously found via scal¬ 
ing and numerical analysis in earlier studies [H13 [ini |TT] 
(a is the vortex core size). Thus, a 2D vortex exhibits 
a breakdown of a linear response to an external force 
/, with a truly nonlinear velocity-force characteristics 
v{f) //|ln/|. 

The “comet tail” texture of a moving vortex leads to 
an elastic energy that diverges linearly with system size 


Ey « ttK {L/^y + \n^y/a), (3) 


with the usual logarithm cut off by the length ^ l/n, 
that diverges with a vanishing velocity. The interac¬ 
tion between two moving vortices strongly depends on 
their velocities and orientation relative to the separation 
vector, r. With the eye on the problems of a vortex- 
antivortex annihilation and nucleation by an imposed 
strain, I find the interaction potential Lt|j _v(^)i for a 
vortex and antivortex moving toward each other, v || r 
(Fig§: 

f4,-v(^) = 27rA: [c-sinh"^(^„/r)] , (4) 

« 27tK I a<^^y <^r, 

\ Inr/o, a«Cr<CCt;j 

and a potential (r) for a pair co-moving with velocity 
V T r (Fig[^: 


Uv.vir) = 2 ttK c-sinh ^{^y/r) + + 1 , (5) 

« 2-kK ■ 


c + r/^. 
in r/a. 


y, a<^^y <^r, 

a r ^y, 


with c = ln^^,/a. Thus, in the annihilation configura¬ 
tion (Figj^, vortex attraction for separation beyond 
is suppressed by the motion. Conversely and even more 


dramatically, I predict that vortex pair motion in the 
transverse configuration (Fi g§ leads to a linear confine¬ 
ment on long scales. 

The above velocity-dependent vortex mobility and in¬ 
teraction qualitatively modify the equation of motion 
for the vortex-antivortex separation. This leads to a 
late-time slowed annihilation dynamics that may be an 
important ingredient in the anomalies observed in the 
experiments[T] . 

Analysis: With the above motivation in mind, I now 
turn to the analysis of moving vortices in a 2D over¬ 
damped xy-model 

jdtO = KV‘^9, with V X V0 = 27ri5(r — r„(t))z, (6) 

searching for a vortex solution d{r,t), that for simplic¬ 
ity I take to be moving at constant velocity defined by 
ry{t) = vt. Despite ignoring a number of ingredients[12j. 
I expect it to be a core description of many systems where 
damping is dominant. 

To this end, I take the solution to be 0{r,t) = 6*„(r — 
vt) -I- 9s{r,t), where 9y{r) = ip = arctan( 2 //x) is the az¬ 
imuthal polar angle that is the standard purely transverse 
solution of the static problem (7 = 0 ), that enforces a 
moving unit of vorticity. The 9s (r, t) part is a nonsin¬ 
gular, single-valued function (with a purely longitudinal, 
curl-free gradient) determined by the requirement that 
9{r,t) satisfies the equation of motion Q. Thus 9s{r,t) 
describes the distortion of a moving vortex about the sta¬ 
tionary form 9(r) = (/), with its spatial Fourier transform 
satisfying 

jdt9sik) + Kk‘^9s{k) = 7 v V0,(k)e-*‘^ "‘. (7) 


The exact solution is easily found either directly for 
9s{r,t) or by first Galilean-transforming to the moving 
vortex frame r' = r — vt, dt — / 9* -I- v • Vr', where the 
distortion is 0'(r',t) = 0s (r' -|- vt,t). 

For a vortex that has been moving forever the Fourier 
transform of the steady state distortion (vanishing for 
= 0) is given by 

, ^ 7v ■ V0t,(k) -2Triky ■ 2 X k 

> Kk^-i-fvk k^k'^-ikyk)^ ^ ’ 


where k„ = jv/K. This leads to the “elastic” energy 
spectrum, |V0(k)p = 47r^(fc^ -|- ky)/[{k.y ■ k)^ -|- k^], that, 
on length scales beyond ^y {k ^ ky) is highly anisotropic, 
akin to that of a smectic liquid crystal. On shorter length 
scales it reduces to that of an isotropic stationary (undis¬ 
torted) vortex[3]. 

In real space the steady-state distortion for a 27r-vortex 
moving along the x-axis, in the vortex frame is given by 
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where x = x/^y, y = y/^v Evaluating above integrals 
numerically and adding the singular part of the vortex, 
9y(v) = (/?, gives the real-space vector fields illustrated in 
Figs. 0 

A transient-state field of a vortex that has been mov¬ 
ing for time t (particularly relevant for the annihilation 
problem) can also be computed exactly and is given by 


Vortex energy: It is of interest to calculate the elastic 
energy Ey = ^ J stored in a moving vortex. In 

steady-state, using (|^ I find: 

Ey = ttK (^\/{L/iyY + 1 -I- ln(^„/a) - sinh“^(^„/L)^ , 


e'^{v,t) = / vVd„(k) 


I — e 






kP' — fkn • k 


• ( 10 ) 


Its key generic features are controlled by three length 
scales ^_L = i/ATt/y, = vt. At time t > t^ = ^y/v, 
such that ^ ^_L ^ 11 , one can see from the solution 

( [To| ) that on scales shorter than an anisotropic domain 
C|| X C-L, the solution reduces to the “comet-tail” steady- 
state one, (§ (Figs. 0 0. On longer scales the vortex 
distortion reduces to d'(r,t) « vt ■ V0y(r'), which when 
combined with the singular part gives 


9'(r',t) « 0y(r') -tvt ■ V0y(r') « d.„(r' -|-vot) = (p. (II) 

Thus on scales outside of the vt x ^Ktj^ domain the 
vortex field reduces to that of an undistorted stationary 
vortex d(r, t) = p at its initial, t = 0 position (see Fig 0 . 
This is a diffusive vortex analog of a “causal horizon” be¬ 
yond which the distortion associated with a moving vor¬ 
tex had not had sufficient time to propagate out. Other 
results (e.g., a vanishing vortex mobility, vortex energy 
and interaction between moving vortices) follow directly 
from the above moving vortex solution. 

Vortex mobility: In the steady-state the power input 
by the external force E to drive the vortex at velocity 
V is balanced by the rotational power dissipated, Pyot = 
J^{dt0){K\/'^0) = f^j(dt0)'^ = J^{d^0)'^, gives the 

vortex drag coefficient, jy = (inverse mobility) [21 

unmn: 
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ylnu. 


Thus, at finite velocity, a previously noted divergence 
with system size L or vortex separation r 1311 [7] is cutoff 
by the velocity-length ~ 1/v, thereby displaying a 
nonlinear velocity-force characteristics i.e., an absence of 
linear response down to a vanishing force. 


that for vanishing velocity, L ^y reduces to In L/a 
of a stationary vortex, but for a rapidly moving vortex, 
L ^ ^y gives the energy @ , that diverges linearly with L 



FIG. 4: Vortex-antivortex separation r{t) as a fnnction of 
time, a solntion of ( |15[ ) (solid), is significantly slowed down 
compared to the “naive” dynamics 'yvdr/dt = —K/r (dashed). 


and with the standard logarithm cut off by the velocity- 
length ^y. This later result is due to the confinement 
of the elastic distortion (that in a stationary vortex is 
uniformly azimuthally distributed) to a comet-tail wake 
of a moving vortex. 


Vortex interaction: To further characterize the nature 
of moving vortices I study vortex-antivortex interaction, 
that strongly depends on their velocities and orientation 
relative to the initial separation vector, ry = r+ — r_. 


Motivated by the vortex-pair annihilation dynamics, I 
first compute the energy E^-^{r\\) = y f 
of a vortex-antivortex pair moving toward each other 
with velocity ±v = iuf along the separation vec¬ 
tor ry. In steady-state the solution is given by 
0v,-v(r,t) = 0 +(r - r+ - vt)-l-0+(r - r+ - vt) -|- 6 > 7 (r- 
r_-fvt)-|- 0 “(r—r_-fvt), with singular (v) and smooth (s) 
components for vortex (at r+) and antivortex (at r_), re¬ 
spectively. The corresponding elastic energy i?v,_v (f II) = 

+ (k^k! - (ky ■ k)2) 

+ (fcCk'„l) 0 ] > is given by 


£^v.-v(r’y) 


2TrK 


-—h In — — sinh ^ — 
Av a ry 


\ + ln(ry/a), a < ry < 


(13) 


where ry (t) = r_|_ — r 


2vt and above I evaluated the asymptotic ry dependence using an approximate 
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hard cutoff ^^/ry on low q. Even for coinciding vortex- 
antivortex positions a linear in system size contribution 
L/^v remains due to elastic energy associated with the 
comet tail of each moving vortex (see Fig.([^). Subtract¬ 
ing this constant self-energy piece I obtain the vortex- 
antivortex interaction, bfH _y(r||) advertised in (j^, that 
is qualitatively weaker and shorter range, falling off as 
l/ry at large separations, ry ^ 

Before moving on, I stress that a full vortex annihi¬ 
lation problem is far richer, requiring analysis of a full 
transient dynamics as vortices accelerate from rest, with 
their velocity-length ^v(t) evolving nontrivially and tails 
limited by the “causal horizon”, growing with t from be¬ 


low to beyond their separation, ry(t). Consequently, the 
nature of the interaction _v(ry,5„) is nontrivially ve¬ 
locity dependent. I analyze the associated dynamics of 
r{t) below. 

Another contrasting geometry of interest is that of a 
vortex-antivortex pair co-moving (see Fig[^ with velocity 
V perpendicular their separation vector r= r^. — r_. In 
steady-state, the solution 0v.v(r,t) = 0+(r — r+ — vt) -|- 
0+(r — r+ — vt) -I- 07(r — r_ — vt) -|- 0“(r — r_ — vt) 
leads to the elastic energy Ev,v(?'_l) = y f = 

^ I ^"fe fc4+(k5k)2 (1 - given by 


A'v,v(?'_l) ~ 27tK 


In — — sinh ^ — 
a r_L 




2ttK 


r±/^v+ ln(C«/a) 
ln(r_L/a), 


a < < r_L, 

a < r_L < 


(14) 


eval uated in the same hard cutoff approximation as in 
(13), and giving ^7.;^.^,(r_L) advertised in (|6|. This is a 


striking result as it predicts for r_L > a linear con¬ 
finement of a moving vortex-antivortex pair, replacing 
logarithmic potential for a stationary pair. As is clear 
from Fig. 0 this elastic energy is associated with the rj_ 
length of the non-overlapping parts of the “comet” tails, 
the rest, beyond r± parts canceling between co-moving 
vortex and antivortex. 

Vortex-antivortex annihilation dynamics, approxi¬ 
mately described (neglecting [T2] transients in ([^) by 
'jydr/dt = — 


dU-v^v{r,v) 


2-kK 


dr 


fln(|f|/2) = 


1 1 


f ^f2f2 _|_ 2 ^ 


(15) 


is significantly enriched|12j by the velocity-dependent 
mobility 0 and interaction 0 , as compared to the naive 
dynamics ^dr/dt = —K/r, that predicts a vortex sepa¬ 
ration r{t) = initially separated by tq, 

annihilating in time to = Above f and 

t are respectively measured in the microscopic units of 
a and ta = a‘^^j(2K). Equation (15) predicts in units 

of Va = a/tg that ^ f ^1 -|- 4/ \\^(v/2) — 1 j 


(rather than rv = const, of the naive dynamics) and 
can be solved numerically, with the result illustrated in 
Fig .(|^. It shows a significant modification and slowing 
of the dynamics by the effects studied here. 

Beyond the transient time the enriched dynamics 
is expected only in the regime of large separation and 
high velocity rv ^ ava = AT/y, corresponding to r >> 
Using AT/y « 10“^ cm^/sec and v = 1/rm/sec I 
estimate « 1mm and Va ~ Imm/sec for a » 1/im, 
a limited regime of current experiment’s[T] applicability. 


Also, above prediction for the product rv decreasing with 
r is inconsistent with measurements [1]. Thus, I conclude 
that in current vortex annihilation experiments, the high 
velocity effects studied here are not sufficient to account 
for the observed anomalies[T] and other effects[Tl] may 
need to be considered. Further systematic experiments 
on moving vortices would be highly desirable to sort out 
various contributions. 

1 also leave the extension of the present London 
limit analysis to a superfluid, beyond a linearized xy- 
model treatment m, incorporating the full Galilean 
invariance m for a future study. 
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